Regional map

We plot a variable of interest from CMAP in a range of latitude and longitude at a given point of time.

Specify following parameters to retrieve required data from CMAP.

• Table Name (table.name)
• Variables to be extracted (sel.var)
• Range variable includes (latitude, longitude, depth, time)

plot_regMap function download the data, and output a list which contain a) regional map plot object: plot_ly/ggplot with “x axis” as “longitude” and “y axis” as “lattitude”; b) corresponding data table.

Example I

con <- dbConnect(odbc::odbc(), DSN="CMAP-MSSQL",UID="ArmLab",PWD="ArmLab2018")

# Inpit variable:
table.name = 'tblsst_AVHRR_OI_NRT'   
sel.var = 'sst'                 
range.var <- list()
range.var$lat <- c(10,70)
range.var$lon <- c(-180,-80)
range.var$time <- c('2016-04-30', '2016-04-30')


## Data retrieval
tbl.subset <- getTableData(con, table.name, sel.var, range.var, order.var=  c('lat','lon'))
head(tbl.subset)

## # A tibble: 6 x 4
##     lat   lon time         sst
##   <dbl> <dbl> <date>     <dbl>
## 1  10.1 -80.1 2016-04-30  29.6
## 2  10.1 -80.4 2016-04-30  29.8
## 3  10.1 -80.6 2016-04-30  29.9
## 4  10.1 -80.9 2016-04-30  29.9
## 5  10.1 -81.1 2016-04-30  29.8
## 6  10.1 -81.4 2016-04-30  29.8

## Plot - regional map 
# out <- plot_regMap(con, table.name,sel.var,range.var, type = 'ggplot')
out <- plot_regMap(con, table.name,sel.var,range.var, type = 'plotly')
out$plot

dbDisconnect(con)

Example II

con <- dbConnect(odbc::odbc(), DSN="CMAP-MSSQL",UID="ArmLab",PWD="ArmLab2018")



table.name <- "tblPisces_NRT"   
sel.var <-  "Fe"               
range.var <- list()
range.var$lat <- c(10,70)
range.var$lon <- c(-180,-80)
range.var$depth <- c(0,0.5)
range.var$time <- c('2016-04-30', '2016-04-30')


## Subset selection:
tbl.subset <- getTableData(con, table.name, sel.var, range.var, order.var=  c('lat','lon'))
head(tbl.subset)

## # A tibble: 6 x 5
##     lat   lon depth time             Fe
##   <dbl> <dbl> <dbl> <date>        <dbl>
## 1    10 -80   0.494 2016-04-30 0.000725
## 2    10 -80.5 0.494 2016-04-30 0.000705
## 3    10 -81   0.494 2016-04-30 0.000658
## 4    10 -81.5 0.494 2016-04-30 0.000635
## 5    10 -82   0.494 2016-04-30 0.000640
## 6    10 -82.5 0.494 2016-04-30 0.000603

## Plot - regional map 
# out <- plot_regMap(con, table.name,sel.var,range.var, type = 'ggplot')
out <- plot_regMap(con, table.name,sel.var,range.var, type = 'plotly')
out$plot

dbDisconnect(con)